Current induced transverse spin-wave instability in thin ferromagnets: beyond linear 

stability analysis 
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A sufficiently large unpolarized current can cause a spin-wave instability in thin nanomagnets 
with asymmetric contacts. The dynamics beyond the instability is understood in the perturbative 
regime of small spin-wave amplitudes, as well as by numerically solving a discretized model. In the 
absence of an applied magnetic field, our numerical simulations reveal a hierarchy of instabilities, 
leading to chaotic magnetization dynamics for the largest current densities we consider. 
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I. INTRODUCTION 

Almost a decade ago Slonczewskii and Bergeip^ pre- 
dicted that when a spin-polarized current is passed 
through a ferromagnet it transfers the transverse compo- 
nent of its spin angular momentum to the ferromagnet. 
The experimental verification of the theoretical predic- 
tions followed within a few years i^i^i^iSii Since then, the 
so-called 'spin-transfer effect' has been observed in a large 
number of different experiments. 

In most experiments, the spin-transfer torque is stud- 
ied in a ferromagnet -normal-metal-ferromagnet tri-layer 
structure where a thick ferromagnet first polarizes the 
current which then exerts a spin-transfer torque on a 
second thinner ferromagnet. At sufficiently large applied 
current densities, the spin-transfer torque then may al- 
ter the magnetization direction of the thin magnet. The 
observation of hysteretic magnetic switching for one cur- 
rent direction only was seen as a hallmark of the spin- 
torque effect^ and excluded an explanation of the ex- 
periments in terms of the Faraday field associated with 
the applied current. (Note that for small system sizes, 
the spin-transfer torque, which scales proportional to the 
current density, dominates over the torque exerted by the 
magnetic field caused by the current flow, which is pro- 
portional to the total current.) Dynamical aspects of 
the magnetic switching process were addressed in recent 
experimentsi2i2ii2iii 

Over the past few years there has been much theoreti- 
cal interest in understanding the spin-transfer torque and 
its consequences for hybrid ferromagnet-normal-metal 
devices. The connection between spin currents or spin ac- 
cumulation in the normal metal spacer layer and the s pin 
torque can be considered understoodiSiiii^ (see Ref. [l^ 
for a recent review). Most calculations of the response of 
the magnetization to the spin-transfer torque have been 
done in the so-called 'macrospin approximation', assum- 
ing that the ferromagnets remain single domains dur- 
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ing spin-transfer induced switching events i^S'il'iSiiSiSfliSi 

They have addressed the precise nature of the magnetic 
switching process, the possibility of limit cycles, and the 
temperature dependence of the spin-transfer torque. In 
addition, full micromagnetic simulations have been done 
by several groups ^SiSiSiiSi e.g., to examine the effect of 
the Ampere field on the hysteretic switching or the break- 
down of the macrospin model into quasi-chaotic dynam- 
ics at very high current densities. While the micromag- 
netic simulations are a significant improvement on the 
macrospin approximation when it comes accounting for 
spatial non- uniformities in the switching process, the ex- 
isting simulations derive the spin-transfer torque from an 
externally fixed spin current, which is a poor description 
of the experimental geometries in which the spin cur- 
rents are determined as an intricate combination of spin 
polarizations caused by all ferromagnetic elements in the 
device-iSiiS-S&Si 

In a recent work, two of the authors showed that a suf- 
ficiently large but unpolarized electrical current flowing 
perpendicular to a single thin ferromagnetic layer can ex- 
cite spin waves in the ferromagnet f2S These spin waves 
have wavevector perpendicular to the direction of cur- 
rent flow. The key mechanism behind the transverse spin 
wave instability is electron diffusion in the normal-metal 
contacts perpendicular to the direction of current flow, 
see Fig. ^ Electrons backscattercd from the ferromag- 
net are spin polarized, the polarization direction being 
antiparallel to the direction of the magnetization at the 
location where they were reflected from the ferromagnet. 
When these electrons reach the ferromagnetic layer a sec- 
ond time, they typically do so at a different point at the 
normal-metal-ferromagnet interface. In the presence of 
a spin wave, the magnetization direction of the ferromag- 
net will be different at that point, and these electrons will 
transfer the perpendicular component of their spin to the 
ferromagnet, thus exerting a spin-transfer torque. The 
sign of this torque is to enhance the spin-wave amplitude. 
A similar argument can be made for electrons transmit- 
ted through the ferromagnet, but their torques tend to 
suppress the spin-wave amplitude. Typically, source and 
drain contacts are asymmetric, and a net spin-transfer 
torque is exerted on the ferromagnet. This torque leads 
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FIG. 1: Through the spin-transfer torque, an unpolarized 
electrical current flowing perpendicular to a thin ferromag- 
netic layer can enhance or suppress spin waves. Electrons 
backscattered from the ferromagnet at point 1 have their spin 
predominantly polarized antiparallel to magnetization direc- 
tion m(l). These electrons exert a torque on the ferromag- 
net's magnetization m(2) if they reach the ferromagnet a sec- 
ond time at point 2, the direction of the torque being to en- 
hance an existing spinwave [i.e., to increase any pre-existing 
difference between m(l) and m(2)]. When electrons transmit- 
ted through the ferromagnet reach the ferromagnet a second 
time at point 3, they exert a torque that suppresses an exist- 
ing spinwave. If source and drain contacts are not symmetric, 
there is a net torque on the ferromagnet, which enhances or 
suppresses the spin wave, depending on current direction. 



to a spin wave instability for the current direction in 
which the effect of backscattered electrons dominates, 
and not for the other current direction. Experiments on 
nanopillars a with single ferromagnetic layer have verified 
the theoretical predictions finding spin-wave instabilities 
for one direction of the current and for asymmetric junc- 
tions onlyiSS 

For a quantitative theory of this transverse spin- 
wave instability, an approach that combines a full self- 
consistent determination of the spin- transfer torque and, 
at the same time, goes beyond the macrospin approxima- 
tion is essential«2£ Indeed, the macrospin approximation 
does not allow for non-uniform spin waves in the ferro- 
magnet, and, whereas an externally imposed spin trans- 
fer torque would predict a similar instability, a non-self- 
consistent theory would be quantitatively incorrect {e.g. 
predict the wrong wavelength for the spin wave) because 
it neglects the coupling between the spin current and the 
spin waves in the ferromagnet. 

The possibility of current-induced non-uniform modes 
in heterostructures has become of recent interest 
in the field, both for single-layer and multilayer 
structuresi^°i^-^i^^i^^i^^i^^ In particular, Ji, Chien, and 
Stilea^S reported experimental and theoretical evidence 
that for large ferromagnet thickness, ferromagnet- 
normal-metal junctions are unstable to the generation of 
non-uniform magnetization modes, but in this case, these 
are longitudinal modes (see also Refs.l^and ls^ . Further, 
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FIG. 2: Spin will accumulate in normal metals on both sides 
of a ferromagnetic layer with uniform magnetization if an un- 
polarized current is passed through the ferromagnet (top left). 
A large-amplitude spinwave in the ferromagnet reduces the 
amount of spin polarization in the normal-metal regions ad- 
jacent to the ferromagnet and lowers the total resistance of 
the device (bottom left). This is shown schematically in the 
circuit diagrams (right). The top two circuit diagrams show 
the resistances seen by majority and minority electrons when 
the magnetization is spatially uniform, the short and long re- 
sistor symbols referring to minority and majority resistances, 
respectively. The ferromagnet with a large-amplitude spin 
wave can be seen as a parallel configuration of ferromagnets 
with opposite magnetization directions. The bottom two cir- 
cuit diagrams show the resistances seen by two spin directions 
in this case. The net resistance is lower in the presence of a 
large-amplitude spin wave. 

Stiles, Xiao, and Zangwill pointed out that transverse 
spinwaves can be excited even in symmetric junctions if 
the spinwave mode is at not uniform in the direction of 
current flow. However, excitation of these modes requires 
a higher currents than the transverse spin-waves consid- 
ered here.'^^ 

Our previous workj^ as well as the other theoretical 
works on this and related spin- wave instabilities;^ii24 was 
a linear stability analysis, sufficient to predict the onset 
of the instability, but not to describe the spin wave ampli- 
tude for current densities larger than the critical current 
density. Knowledge of the spin wave amplitude is nec- 
essary if one wants to study, e.g., how the spin wave 
instability affects the resistance of the normal-metal- 
ferromagnet junction. It is the goal of this present work 
to examine in detail the dynamics of the spin-wave be- 
yond the instability. While we focus on the case of single- 
layers, we expect that, in light of the work of Refs. l32l34 . 
our qualitative findings will carry over to the case of tri- 
layers and heterojunctions. 

Although a quantitative description of how the spin- 
wave instability affects the resistance of the normal- 
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metal-ferromagnet junction will be postponed to the 
next two two sections, the sign of the effect can be deter- 
mined using simple considerations. Once the current den- 
sity has exceeded the critical current density for the spin 
wave excitation and a spin wave has been established, 
the fact that the magnetization is no longer uniform re- 
duces the amount of spin accumulation in the normal 
metal contacts adjacent to the ferromagnet. A reduction 
of the spin accumulation in the normal metal contacts 
causes a reduction of the sample's resistance, see Fig. 
121 for a schematic drawing. Indeed, the experiments of 
Ref. observed a small decrease of the resistance of the 
nanopillar upon onset of the spin-wave instability. The 
effect of a purely transverse spinwave instability is oppo- 
site to that of a longitudinal spinwave, which increases 
the resistance of the device. ■^^ The reduction of the spin 
accumulation in the normal-metal spacer also lowers the 
spin-transfer torque, thus providing a mechanism to sat- 
urate the growth of the spin wave amplitude for current 
densities larger than the critical current density. More- 
over, note that a theory of this effect needs to combine 
features of both the micromagnetic approach and the self- 
consistent treatment of the spin-transfer torque. 

In Sec. ini we consider current densities slightly above 
the critical current density. In this regime, a perturba- 
tive treatment in the spin wave amplitude is possible. 
In Sec. mil we then perform a detailed numerical sim- 
ulation of a simplified system that allows us to probe 
current densities much larger than the critical current 
density. Whereas the observed magnetization dynamics 
in the presence of a large magnetic field is rather unsur- 
prising — there is one stable energy minimum, and the 
magnetization precesses around the direction for which 
energy is minimal — , in the absence of an external mag- 
netic field we find a hierarchy of instabilities. For very 
high currents the system shows chaotic behavior with 
measurable Lyapunov exponents. 



II. PERTURBATIVE CALCULATION 

We consider a single ferromagnetic layer, connected 
to source and drain reservoirs, see Fig. 13 Between the 
ferromagnet and the drain reservoir is a normal-metal 
spacer, as is common in nanopillar geometries. There is, 
however, no normal-metal spacer between the ferromag- 
net and the source reservoir. We use coordinates x, y, 
z, where x is the coordinate perpendicular to the layer 
structure and y and z are coordinates in the plane of the 
layers. 

Both the ferromagnet and the spacer layer have a rect- 
angular cross section of dimensions Wy x Wz ■ The ferro- 
magnet has thickness d, which is taken small enough that 
the chemical potential for the conduction electrons and 
the the direction m of the magnetization of the ferromag- 
net do not depend on the longitudinal coordinate x. The 
normal metal spacer has thickness L. Transport through 
the normal metal spacer is diffusive, with conductivity a. 
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FIG. 3: Schematic picture of the normal-metal-ferromagnet- 
normal-metal junction considered in our calculations. The 
ferromagnetic layer (F) is connected to source and drain reser- 
voirs though normal metal spacers (N) . We consider the maxi- 
mally asymmetric case with only one spacer of length L 3> ^sf . 



In the normal metal spacer, the charge and spin de- 
grees of the conduction electrons are described by the 
equations 



V^Aic = 0, 



(1) 



where and fi^ are chemical potentials for the electron 
density and electron spin respectively, — e is the electron 
charge, and Igf is the spin diffusion length in the normal 
metal spacer. Further, is the charge current density 
and a is the conductivity of the normal metal leads. The 
boundary conditions for x — L at the drain reservoir is 



^i,iL)^-eV, At,(L) = 0. 



(2) 



Here the argument L refers to the x coordinate. The y 
and z coordinates are not written explicitly. The second 
boundary is the interface between the normal-metal and 
ferromagnet at x — 0. Since the electron dynamics hap- 
pens on a time scale that is much faster than the rate 
of change of the magnetization direction m, this bound- 
ary condition can be taken treating m in the adiabatic 
approximation 
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[.g+Aic(0)-l-5-m-/i3(0)] 



-^T (5-Mc(0) + .g+m • /x,(0)) m 



+ ^ffi (2/is(0) X m + hm) x m 

+ ^^2(2^^(0) xm + ririi). (3) 

Here g± — {g^ ±gi)/2, where g^ and gi are interface con- 
ductivities for spins aligned parallel and anti-parallel to 
m, and gi -t- ig2 is the complex valued 'mixing interface 
conductivity'. The argument "0" refers to a coordinate in 
the normal metal spacer, just outside the ferromagnetic 
layer. The charge current and the spin current parallel 
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to m are continuous at the interface. In writing down 
Eq. (O we assumed that the two ferromagnet-normal- 
metal interfaces are identical, so that the potentials 
and m • fi^ drop equaUy over both interfaces of the ferro- 
magnct; the component of /Xg perpendicular to m is zero 
in the fcrromagnct. (It is the non-conservation of spin 
current perpendicular to m that gives rise to the spin 
transfer torque.) For Co/Cu and Fe/ Cr inte rfaces, these 
conductivities are tabulated, see Refs. ll4l37l Typical val- 
ues are in the range 92 gi ^ g± ^ 10^"* n~^ni~'^. For 
any interface, one has the constraint gi > g+ > 5-.^^ 

We are interested in the situation in which the magne- 
tization is allowed to vary in the direction perpendicular 
to the current flow. In this case a large enough current 
may cause spin-wave excitations perpendicular to the di- 
rection of current flow.^^ To simplify the notation, we 
take the limit L 3> ^sf • The spin and charge chemical po- 
tentials in the normal-metal spacer then have the general 
solution 



,(r) = ^e*««^+*«^^ac(q)e"^'^« 



elx 



Ms(i-) = ^e*««^+*«^^e-('2+«'+'^^')'^'^a,(q), (4) 



where q = {0,qy,qz)'^ is a wavevector in the y-z plane. 
The components qy and qz take values qy = nny/Wy, 
qz — nriz/Wz with integers Uy and Uz. The Fourier ex- 
pansion coefficients ac(q) and as(q) are real and satisfy 



ac(q) ac(-q), a^{q) = as(-q) 



We further define the quantities 



(5) 



^2n1/2 



Gc(q) - {<j/2){ql 
GM) = {<^/2){ql+ql + Q')'/', 



(6) 
(7) 



which have the same dimension as the interface conduc- 
tivities g±, gi, and 32 • With these definitions, the bound- 
ary condition Q at the normal-metal-ferromagnet inter- 
face becomes 



el 



WyWz 



2E' 



,iqyy+iq^z 



[Gc(q)ac(q) + g+a^iq) + g-a^iq) ■ m] , 



[2Gs(q)as(q) + (g_ac(q)m -Hg+as(q) • m)m - 25i(as(q) x m) x m - 2g2as(q) x m) 



hgiui x m — hg2Yn.. 



(8) 



Although Eq. I© gives a set of linear equations for 
the expansion coefhcients ac(q) and as(q), a solution in 
closed form is not possible for arbitrary magnetization 
m(2/, z). Instead, we expand around the uniform equilib- 
rium direction. Hereto we introduce a second coordinate 
system with axes labeled 1,2, and 3, such that m points 
along the unit vector 63 in the absence of an applied 
current, and write 

m = miei H- 771262 -f- (1 - - 7712)^/^63. (9) 

We then perform a Fourier transform, similar to Eq. (0J 

m,(y,z) = ^m,(q)e'««''+^«^^ ^' = 1,2, (10) 
q 

where TOj(q) = mj[—q). Finally, expanding in powers of 
mi and TO2, we have solved the spin and charge chemical 
potentials to third order in mi and m2, which parame- 
terize the deviations from equilibrium. 

In order to complete the calculation, we need to cal- 
culate the rate of change of the magnetization direction 



m in the presence of the current /. Hereto we use the 
Landau-Lifshitz-Gilbert equation^^Si^ 

m = am X m + Tcx + Tan+Tno, (11) 

where a is the Gilbert damping coefficient, Tcx is the 
torque arising from exchange. Tan is the torque from the 
combined effect of magnetic anisotropy and an applied 
magnetic field, and Tno represents the current-induced 
spin-transfer torque. The latter reads^ 

T-„e = -^(js(O) -js(-rf)) X m) X m 

= [51 (/^s X m + /irh) X m 

-H52(Ats X m + ftm)]. (12) 

Here the spin current js{—d) is taken in the source reser- 
voir, M is the magnetization per unit volume and 7 is the 
gyromagnetic ratio. Note that the terms proportional to 
the time derivative m have contributions from two inter- 
faces while the contribution to the torque from the spin 
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chemical potential has a contribution from the a; = in- 
terface only. (All potentials are zero in the source reser- 
voir.) The exchange torque Tox is 



JjMV^m X m, 



(13) 



where J is the exchange constant. To linear order in mi 
and m2, the anisotropy torque Tan can be written 



— (fcimiGi 



fc2m2e2) X m, 



(14) 



where ki and ^2 describe the combined effect of magnetic 
anisotropy and an applied magnetic field. If anisotropy 
dominates over the effect of a magnetic field, higher-order 
terms in an expansion in powers of mi and m2 will be 
highly sample specific. Although this case can be dealt 
with using the methods presented below, the result of the 
calculation has little predictive value if those coefficients 
are not known independently. Therefore, we focus on 
the opposite limit that the anisotropy term in Eq. H14|l 
is dominated by magnetic field. Then higher-order terms 
in an expansion in powers of mi and m2 are related to 
the first-order terms, and one has 



Tan = (fc7/Af )e3 X m. 



(15) 



where we wrote ki = k2 = k. For future reference, we 
combine the material constants J and 2k — ki + k2 into 
the combinations 



(16) 



which have the dimension of inverse length and current 
density, respectively. 



We now proceed to report the result of our calculation. 
The lowest order result, indicated by a superscript "(0)" , 
is 



„(0)(q) ^ ej(5+ + 2Gs(0)) 



a(°)(q) 



e3(5q,o. (17) 



5m (0) 

Here j = I/WyWz is the current density and2& 
. ^ (g+ + 2G,(q))(ff++2Ge(q)) 

5m q = g-- 18) 

9- 

Writing ^c(i) = -eV = -e{L/aWyWz + R)I, we con- 
clude that the resistance R of the ferromagnetic layer is 



R 



WyW, 2g+a/kf+gl 



2 ' 

9- 



(19) 



For the zeroth-order solution, the spin potential /ig is 
collinear with m throughout the sample. Hence, to that 
order there is no current-induced torque. This is different 
when small deviations from the situation m = 63 are 
taken into account to first order. One finds that the first- 
order corrections ai^^ (q) and a^^^ (q) are zero. In order to 
represent the first-order contributions to the transverse 
spin potentials Ogi and as2, we use spinor notation, = 
(asi,as2)"'" and m = (mi,m2)'^. Then, defining 



^(q) = (.9i + G.(q))^ 



52. 



(20) 



we find 



'5m(0) 

h 



m(q) 



ej(Gs(q)-Gs(0)) 
5m(0)i?(q) 



[(51 + Gs(q))m(q) + ig2cr2m{q)] 



[52Gs(q)m(q) + {gf + giGs{q) + gi)ia2m{q) 



(21) 



where (T2 is the second Pauli matrix. Note that the first 
term on the right hand side is the response to a uni- 
form rotation of the magnetization, while the second and 
third terms give the response to a non-uniform and time- 
dependent magnetization. 

The potentials are substituted into Eq. 1)12(1 to find 
the current-induced torque, and then into the Landau- 
Lifshitz- Gilbert equation (|ll|l to find the rate of change of 
the magnetization. The current-induced torque has con- 
tributions proportional to the time derivative m, which 
lead to a renormalization of the Gilbert damping param- 



eter a and the the gyromagnetic ratio 7. The renor- 
malized Gilbert damping parameter a and gyromagnetic 
ratio 7 = 7//? depend on the transverse wavevector q 
and read 



1 



2Mde^ 



1 



Diq) 



2Mde^ 



1 



Diq) 



(22) 



In the macrospin limit q —> 0, these modifications coin- 
cide with the renormalized values originally reported in 
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Ref.|27l 

Again using two-component spinor notation, the com- 
plete Landau-Lifshitz-Gilbert equation then becomes 



012 + ia2a)m{q) = A(q)m(q), 



(23) 



with 



^(q) 



1(72 



(1) , ^7Jf(g^ +gf ) 



2eqfM 



(T3- 



7(^1 - h) 



and 

^(1) 



II 

■f(q) 



2M 



9l 



(24) 



Mde^ 5m(0)i?(q) 
Mde^ .gm(0)i?(q) 



[G,(q)-G,(0)], 
[Gs(q)-Gs(0)]. (25) 



In the absence of a current, any spatial modulation 
of the magnetization is damped. However, a sufficiently 
large positive current / can overcome the damping, and 
cause a spatial modulation of m to grow in time, rather 
than decay. (A positive current / corresponds to elec- 
tron flow in the negative x direction.) The instability 
condition is easily obtained from Eq. H23|l 



<5(q) 



2eqfM 



+ r«(q). (26) 



We can analyze this result in different limits. For a ferro- 
magnetic layer with sufhciently small transverse dimen- 
sions, Wy,Wz < (Isf/qf)^^^ if IsiQf > 1, the instability 
happens at wavevector q = {n/Wy)y or q = {tt/Wz)^, 
whichever is smallest, and the critical current follows di- 
rectly from Eq. (|26l) . For wider layers, the critical cur- 
rent density jc and critical wavevector qc are found as 
the current-density wavevector pair for which the onset 
of the instability condition happens at the lowest current 
density. 

This condition can be simplified in the limit of a very 
thin ferromagnetic layer, d — > 0, neglecting terms propor- 
tional to g2 (which is numerically smaller than gi), and 
for wavenumbcrs g <C Qf. We then find that the critical 
current follows from minimizing the relation 



ic(q) 



'75m(0)jf 



q + qf 



2Mqfe^ I - [I + qH^J- 



1/2- 



In the limit ^sf 3> l/qi, this gives^S 



^^7gm(0)gfjf 

2Me2 



(27) 



(28) 



(The result for jc was reported incorrectly in Ref. I23. 
Note that the condition qc <C gf , which was used to derive 
Eq. 123 is consistent with Eq. J^Hl) if ki > l/ft-) Note 
that qi increases with an applied magnetic field, so that 
this limit becomes relevant even for the case of a normal 



metal with strong spin relaxation if the magnetic field 
is large enough. In the limit Igf ^ 1/qi of strong spin 
relaxation and weak anisotropy, one has 



= (4/3)i/4(gfA.t)i/2, Jc ^ 



Mqil^^fe^ 



(29) 



At the critical current density, the trajectory of the 
magnetization is a simple ellipse (circle in the case of 
large magnetic fields) . The ellipse is described by the co- 
ordinate transformation mi — r(cos 6 cos (p+rj sin 9 sin (p), 
and 7712 = r {sin 9 cos (p — 77 cos sin 0). The solution of 
the magnetization dynamics at the critical current then 
gives — ujQt and r constant, where ujq — lu'^ — uj'^, 
77 ~ {uj+ — uj- ) I Wo and 

_i _ 2Mce^qiCOs(29) 



7(9f + ql)k 



2Mce^qfg2Gsici) sin{29) 



2Mch 



2 
92 



-5iG.(q)))' 



7(^1 



(30) 



(31) 



and c, 9 are obtained from a = csm29, $ = ccos29. For 
the case of a large applied magnetic field, ki ^ k2 = k, 
and neglecting g2, we have i] = 1 and 



wo = ^7if('?c+'?f)/(2egtM). 



(32) 



Note that, although the applied current has a large ef- 
fect on the stability of the ellipsoidal motion (precession 
is damped for j < jc and unstable for j > jc), its effect on 
the precession frequency is small. To a good approxima- 
tion, the precession frequency equals the ferromagnetic 
resonance frequency in the absence of a current. 

Whereas the first-order calculation allows one to find 
the current density at which the spin-wave instability sets 
in and the angular form of the low-amplitude excitations, 
it does not provide information about the magnitude of 
the spinwave oscillation for j > jc , or about the effect of 
the spinwave oscillation on the resistance of the ferromag- 
netic layer. This information can only be obtained from 
the analysis of the magnetization dynamics beyond first- 
order in the amplitude. Such a program proceeds along 
the same lines as the first-order calculation shown above: 
Calculation of the potentials for charge and spin in the 
presence of a non-uniform and time-dependent magneti- 
zation, followed by a calculation of the current-induced 
torque and the rate of change of the magnetization. We 
have carried out this program to third order in 7771 and 
7772 , and list some of our general results in the appendix. 
However, as this calculation involves higher-order con- 
tributions to the anisotropy torque Xan, for which the 
expansion constants are unknown, we find that this cal- 
culation has little predictive value. Instead, we focus on 
the limit in which all magnetic anisotropy arises from an 
applied magnetic field. In this limit. Tan is known, cf. 
Eq. H15|) . and a theoretical analysis is useful. 
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An important simplification is that the higher-order 
analysis is necessary for the Fom'ier components mi(qc) 
and m2(qc) at the critical wavevector only. Hence, we 
need to consider only a single Fourier component in our 
considerations below. Solving for the leading (second or- 
der) correction to the charge potential, we find an ex- 
pression that depends on the magnetization amplitude, 
to second order in mi and rn2, and on the time deriva- 
tives. Only first-order time-derivatives appear, which can 
be eliminated using the Landau-Lifshitz-Gilbert equation 
For the case of a large apphed magnetic field, the 
magnetization precession is circular, and one has 

mi(qc)m2(qc) - ™2(qc)™i(qc) = t^o»'(qc)^, (33) 

where we abbreviated 

r(qc)2 = mi(qc)2 + m2{qcf . (34) 

The precession frequency ujq given by Eq. H3U|) above. We 
then find 

_ 2(G,(0) - G,(qe))r(qe)^ 
' ^ ^(qc)5m(0)2 

X [u;o5m(0) (i)(qc) - G,(qc)(5i + G.(qc))) 

- 2eje (D(qc) + (Gs(0) - Gs(qc))(.9i + Gs(qc)))] 

(35) 

Solving for the leading (third) order torque, we note 
that the third order torque depends not only on the mag- 



netization amplitudes mi{qc) and 7n2(qc), but also on 
their time derivative rhi and 1712 ■ The time derivatives 
appear to first, second, and third order in the expansion. 
The dependence on m'^^^ leads to the same modifications 
to the Gilbert damping and gyromagnetic ratio as for the 
first-order current-induced torque calculated above. The 
dependence on m^^^ is through the 3-component only, 
which can be written as 



,(2) 



(1) 



• (1) 

77127712 . 



(36) 



The first-order time derivatives mf^^^ can be expressed 
in terms of 7711 and m2 using Eq. I|33(l [or, in the gen- 
eral case, using Eq. (|23|l ]. For the anisotropy torque Tan 
we take the contribution from the magnetic field only. 
Hence, 



Tcx(qc) + T-an(qc) 



X (-«cr2)TO(qc). 



yKqc 



(37) 



Thus proceeding, we find that the third-order equation 
for the rate of change of the magnetization direction reads 



(/3I2 +i(T2(5)rh(q)(3) = ^(q)(3)m(q)(3), 



(38) 



with 



J 



^(qc 



U3) = 



2T| ^^(0)l2 - 2ia2Tf\0) + T ^^(2qc)l2 - ifTzT^' (2qe) + 35c^ol2 



2eq{M 



ia2 



(39) 



and 



rf(k) 



%gi [ [9+ + 2Gc(k)] [g+ + 2G3(q)] - 
Mde^ig, + Gq) \ [g+ + 2Gc(k)][.g+ + 2G,(k)] - 



5m (0) 



[Gs(k) - G,(0)] 



gihuf) - 



2ej(.9i + G,(0)) 



5m (0) 



G,(k) - Gs(q) 
.9l + Gs(q) 



hjg2 G,(q)[Gs(k)-G,(q )] f [.9+ + 2Ge(k)][.g+ + 2G,(q)] - g^_ 
Mde^gi+Gs{qW 



' ' , ej(<7i + Gs(0)) 



[.9+ + 2Gc(k)][5+ + 2G.(k)] - gi ["^'^"^ 5m(0) 
[{g+ - 2gi)(2Gc(k) +g+) -ff2][G^(q) _ Gs{0)]ej 
gmmi9+ + 2Gc(k))(.9+ + 2Gs(k)) - gl] 



(40a) 



(40b) 



Solving the differential equation for ?7i, one finds that the precession amplitude for current density j slightly above 
the critical current density jc reads 



5(9? + 9f ) 



eqfMjc p 



2r|f)(0) + r|f)(2qe)-3T|J^^(qe) 



2Tf ^(0) + rf ^(2qc) - hjjtq^J {2eqfM) 



(41) 



The result takes a simpler form in the limit 172 — (since g2 is numerically smaller than gi), d — > 0, and l/lst 
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Qc < q{, 



(j -Jc)(g^ 



2 

9- 



3ci'2gig+ +g+a/ls 



(42) 



Since gi > g+ > g^ we conclude that the r(qc)^ > is 
positive if j > jc, which excludes hysteretic behavior. 

In the same limit we can also calculate the change in 
frequency of the spinwave given by 



^ = 1 



(43) 



Since the prefactor of the second term is much smaller 
than unity, Qc <C qt for the parameter regime of interest, 
we conclude that in the regime of perturbation theory, 
there is hardly any change from the ferromagnetic reso- 
nance frequency. 

Finally, at the onset of the spin-wave instability, the 
resistance of the ferromagnetic layer acquires a small neg- 
ative correction 



R_ 

i?o 



(0) 



«^°^(0)' 



1 - 



2((T/?sf + 351)52 r(qe)2 



gi + g+cr/lsf) ' 



(44) 



(In the second line we took the limits 52 ^ 0, d ^ 0, 
and used l/lgf ^ qc ^ ^f-) This resistance decrease is 
anticipated on physical grounds since the non-uniform 
mode allows for an increased transmission of minority 
elections that diffuse along the transverse direction — 
see Fig. Inland the corresponding discussion in SecQ] 



III. NUMERICAL CALCULATION 

The calculations in the preceding section are valid for 
currents close to the onset of the instability. For currents 
much larger than the critical current, we need to go be- 
yond perturbation theory to obtain the dynamics. Hereto 
we numerically solve for the magnetization dynamics and 
its effect on the resistance of the ferromagnetic layer. 

In our numerical analysis, we assume Wz <^ Wy and 
impose that the magnetization direction ni{y, z) does not 
depend on z. The remaining two-dimensional problem is 
replaced by a finite number of one-dimensional problems 
by substituting the normal-metal spacer and the ferro- 
magnetic layer by N normal metal channels, each at- 
tached to a magnet with magnetization direction m(n), 
n — 1, . . . , N. In order to model a higher-dimensional 
structure, electrons are allowed to diffuse between the 
channels, whereas the N magnets interact via an ex- 
change energy. A schematic drawing of this model is 
shown in Fig. 01 

In this discretized model, the potentials for charge and 



v = o 



x = 



x = L 



FIG. 4: Schematic drawing of the model solved numerically. 
The continuous magnet is replaced by magnets (left) , each 
coupled to a normal- metal wire (right). The wires are coupled 
via transverse diffusion (shown schematically as solid lines); 
the magnets are coupled via the exchange interaction (shown 
schematically as dashed lines). 



spin obey the equations 



2iic{n, x)] ~ 0, 



N 



dxtJ'sin,x) + 77^ [fi^{n+l,x) + fj,,{n-l,x) 



2^l^{n,x)] = 



/x,(n,x) 



/2 



(45) 



Equations for the boundary channels, n = 1 and n = 
N, are obtained by setting /j,c,s(0,a;) = pLc^sO-jX) and 
IJ-c.s{N -I- l,a;) = ^c,x{N,x). The general solution of Eq. 
()45|l is of the form 



^ic{n,x) = 2 ^ ac(0cos[/7r(n+ l/2)/A^], 



elx 



aWyWz 



ti,{n,x) = 2 ^ a,(0cos[;7r(n-M/2)/iV]e~«=(')^, 



i=0 



(46) 



with 



9c(0' - 4{N/Wyfsm^ln/2N)), 

qSf = 1;^^ + 'iiN/Wy f sin' {lTr/2N)) . (47) 

The boundary conditions at a; = (normal-metal- 
ferromagnet interface) are given by Eq. |(2Jl. 

The magnetization dynamics is given by the Landau- 
Lifshitz-Gilbert equation with a discretized ex- 

change torque Tex, 

, , J-iMN'^ , , 
'^ex[n) = — — [m(n -I- 1) -I- m(n - 1)J X m(n). 



VVy 



(48) 
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For the anisotropy torque we consider two different cases: 
The hmit of a large appHed magnetic field, 

Tan = -^63 X m{n), (49) 

as well as the case of no applied field, where we take a 
simple model for the torque arising from magnetocrys- 
talline and shape anisotropy, 

Tan(n) = - — [/cimi(n)ei + fc2"i2(n)e2] X m(n). 

(50) 

The Landau-Lifshitz-Gilbert equation, together with 
the boundary conditions at x = 0, are sufficient to de- 



termine the AN expansion coefficients ac{l) and as{l), 
I = 0, . . . , N — 1, and the time derivative of the magne- 
tization directions m(n), n = 1, . . . , iV. Our numerical 
procedure consists of first expressing m(n) in terms of 
the potential expansion coefficients ac(/) and as{l) using 
the Landau-Lifshitz-Gilbert equation, and then solving 
for the potential expansion coefficients using the bound- 
ary condition at a; = 0. 

For the practical implementation of this scheme, it is 
useful to define 3x3 matrices and TZ = mm^ such 
that for any vector v, v x m = A4v and TZ — I3 — 
M.^ . In 3 X 3 matrix notation, the time derivative of the 
magnetization vector can be expressed in terms of the 
potential coefficients as 



m(n) 



(i'tz + a'^n/f3' - a'M 

X < X[rex + Tan] 



2^7 
Mde' 



[Mgi + <?2l3]Xa,(0 cos Mn + 1/2) /N] \ , 



(51) 



where a' = a + h^^gi/ {Mde^) and (3' = 1 + t^^gij {Mde"^). In turn, the potential coefficients a.s{l) are obtained from 
inverting a 4iV dimensional matrix equation, 



2cos(?7r(n + l/2)/iV) 



1=0 

where we abbreviated 



2crgc(0 + 2.g+ 2m'^.g_ 
2mg_ aqs{l)t3 + ixi 



a,{l) \ ^ 2el 
as(/) J WyW^ V S-m 





X2 



(52) 



Xi{n) = g^n-M{giM + g2t3) 

+ 2Mde^\j^ + P'^) ^^^^ + 52l3)[/3'l3 + c^'^n/P' - a' M]M{giM + 52I3) 

X2(") = ^Mr2(t^?+'jf/ + C."nlP' - a'M]M[TM + Tan(n)]. 



(53) 
(54) 



We have performed numerical simulations for N rang- 
ing between 10 and 20, although all data shown are for 
A'' = 10 and N = 11. We verified that there is no 
qualitative dependence on the parity of N in our sim- 
ulations. A small random torque was added at each time 
step to mimic the effect of a small but finite tempera- 
ture. (The corresponding temperature obtained from the 
fluctuation-dissipation theorem was less then a niK.^^) 

Below we present our results. We first consider the 
case in which the anisotropy torque is dominated by an 
applied magnetic field, taking Eq. H49|) for the anisotropy 
torque Tan- We then consider the case in which there is 
no apphed magnetic field, taking Eq. (|5n|l for Tan- The 
latter case is qualitatively different from the former, as 
it has two stable equilibria for m (m = 63 and m = 
—63), whereas in the presence of a large applied field the 



equilibrium position is at m = 63. 



A. Large applied magnetic field 

For the numerical simulations with a magnetic field, 
we took values for the various parameters as follows: 
thickness d = 0.2 nm. Width Wy = 55 nm, as is 
appropriate for typical nanopillar experiments^ spin- 
diffusion length = 100 nm, a/kf = lO^^n-^m-^, 
gi = 5.5 X 10" n-^m-^, 32 = 0.3 x lO" ^-^m-^, 
5T = .9+ + 9-_ = 4.2 X 1014 fl-^m-^, g^ = g+ - g_ = 
3.3 X lO^** O ^m ^. The interface conductivities are 
taken from numerical calculations for a disordered Cu/Co 
interfaceji^ the conductivity a and the spin relaxation 
length Isf are consistent with those in Cu. We further 
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12 3 4 5 

FIG. 5: Main panel shows the magnetization component 
7713(1) of the first magnet, as a function of applied current. 
The solid line is obtained from the perturbation theory result 
1411 1. while the dashed line is a guide to the eye. In a large 
magnetic field, the motion is circular, an example is shown 
in the inset where j = l.Sjc- 

took a = 0.01, h'jgi/Mde'^ = 0.0138, if = IQi^A/m^, 
qt = lO^^nm"^ (as is appropriate for Co, see Ref. liH : 
the magnetic field corresponding to the values of j{ and 
qt listed above is of a strength comparable to the intrinsic 
anisotropy energy). For these parameters, the width of 
the sample is so small that the spinwave wavenumber q 
is set by the finite sample width, q = n/Wy. 

For current densities below jc, no spinwaves are ex- 
cited. Simulation runs in which the magnetization is 
tilted away from the easy axis 63 show damped pre- 
cession towards the equilibrium magnetization direction 
m = 63. For current densities above jc, a- spin- wave with 
wavenumber q — Tr/Wy is excited. Each magnet n in 
our simulation n = l,...,N shows circular precession 
around the direction of the applied magnetic field, see 
Fig. El inset. The amplitude of the oscillation increases 
with current as predicted by the perturbation theory of 
the preceding section. The 3-component of the magneti- 
zation is a constant of the motion and can be monitored 
to measure the amplitude. Numerical results for for 
the magnet n = 1 are shown in Fig. |S1 as a function of 
current density, together with a comparison of our nu- 
merical results with the perturbative result (|41|l . With a 
large applied field, the magnetization dynamics remains 
regular even for current densities much larger than jc- 
The effect of the spin-wave instability on the resistance 
of the ferromagnetic layer is shown in Fig. |S1 

B. No applied magnetic field 

We have also performed numerical simulations in the 
absence of an applied magnetic field. Hereto, we choose 
Eq. H50|l for the anisotropy torque, and choose fci and 
k2 such that (fci - fc2)/(fci -I- = 0.99. This form of 
the anisotropy is appropriate for thin magnetic layers, in 
which the magnetic anisotropy is predominantly of easy- 
plane type. The magnitude of the anisotropy energy is set 




1 2 .„ 3 4 5 

FIG. 6: Resistance of the ferromagnetic layer, as a function 
of applied current (crosses). The solid line is obtained from 
the perturbation theory result l|44|l. while the dashed line is 
a guide to the eye. 



by the parameters qt and jf , for which we take the same 
values as in the previous subsection. All other parameters 
are also taken the same as in the previous subsection. 

The magnetization dynamics without applied magnetic 
field is much richer than the magnetization dynamics at 
a large magnetic field. The reason is the existence of 
two stable equilibrium directions if no external magnetic 
field is applied (m — 63 and —63). At sufficiently large 
current densities, the current-induced torque drives the 
magnetization direction between these two stable direc- 
tions, leading to a variety of dynamical phases. 

For the numerical parameters chosen in our simulation, 
we observe the following characteristic dynamical modes: 
For current densities jc < j ^ 2jc the instability develops 
with the wavenumber q = n/Wy. Because the magnetic 
anisotropy energy used for the simulation has no rotation 
symmetry around the 3 axis, the magnetization direction 
m(n) of each magnet n = 1, . . . ,N traces out an ellipse, 
rather than a circle. We describe the magnetization mo- 
tion is described using Poincare sections for the polar an- 
gles 9 and (f) for the magnetization. The top right panel 
in Fig. 13 shows traces that are symmetric about (j) = t:, 
which have the functional form for m as predicted by the 
perturbation theory in the preceding section. 

For higher currents with 2jc 3 2.5jc, the refiection 
symmetry about the easy axis is spontaneously broken, 
resulting in asymmetric ellipses (upper inset in Fig. |SJ), 
which for even higher current densities turn into orbits 
around the direction perpendicular to the easy axis (lower 
inset in Fig. ISJ. A three-dimensional rendering of this 
regime is shown in Figure |S1 

For even larger currents there is a transition into non- 
periodic modes that cover a significant part of phase 
space, as shown in Figure |21 Whereas these modes are 
non-ergodic for lower current densities, they eventually 
become ergodic and chaotic at high current densities, 
with Lyapunov exponents increasing with the current 
density j (data not shown). 

In this general case, when the magnetization motion is 
not just simple circular precession, the spin-wave insta- 
bility not only leads to a decrease of the dc resistance of 
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FIG. 7: Typical elliptical trajectory for one of the discrete 
nanomagnets m(n) for weak easy axis and strong easy plane 
anisotropy with jc < j < 2jc (left panel). The upper and 
lower right panels show the corresponding Poincare sections 
for j = 1.2 jc and 1.5jc respectively. This regime agrees with 
the perturbative calculation of Sec. ^ where the lowest en- 
ergy spin-wave mode is excited and increasing the current 
only changes the amplitude of elliptical oscillation. 
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FIG. 8: First manifestations of further dynamical instabilities 
in the range 2jc < j < 2.5 jc- The upper right panel shows a 
Poincare section for j — 2.2jc where the motion is no longer 
symmetric about the easy axis. The lower right panel shows 
the motion for j — 2Ajc where the motion is trapped between 
the ±63 easy axes direction. The left panel shows what this 
motion looks like on the unit sphere. 



the ferromagnetic layer, it also causes a fast oscillation 
of the resistance as shown in the time trace in Figure ITUl 
The right panel in Fig. EH shows the decrease of the dc 
resistance up to j = 2.5jc- (No sufhciently accurate nu- 
merical results were obtained for larger current density.) 
Results for the variation of the resistance amplitude and 
frequency with the applied current density are shown in 
the left panel for current densities up to 4jc. At the pa- 
rameter values considered in our simulation, the onset 
of the non-periodic magnetization variations is accompa- 
nied by a sharp rise in precession frequency and a de- 
crease of the amplitude of the resistance fluctuations. 



IV. CONCLUSION 

We have presented a detailed study of the transverse 
spin-wave instability for a single ferromagnetic layer sub- 
ject to a large current perpendicular to the layer. Our 
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FIG. 9: Poincare sections for the magnetization direction of 
one of the magnets at j = 2.5jc (left) and j — 3.2jc (right). 
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FIG. 10: The upper panel shows the time trace of resis- 
tance where the spin- wave instability causes a decrease in the 
observed resistance. The lower left plot shows how the ampli- 
tude and period of the resistance oscillation change with the 
driving current, while the lower right panel shows the decrease 
of dc resistance. 



calculations have been in the small-amplitude regime, 
where perturbation theory can be used, and in the large- 
amplitude regime, where the magnetization dynamics 
can be solved numerically. 

The two main signatures of the spin-wave instability 
are (1) existence of the instability for one current direc- 
tion only, and (2) a small reduction in the dc resistance 
of the ferromagnetic layer. The resistance decrease arises 
because the existence of a spin wave with large ampli- 
tude lowers the spin accumulation in the normal metal 
adjacent to the ferromagnet. A lower spin accumulation 
corresponds to a lower resistance (just as a high spin 
accumulation state of the antiparallel configuration in 
the standard current-perpendicular-to-plane giant mag- 
netoresistance geometry gives a high resistance state). 
Both features have been seen in a recent experiment. 

An important question for a dynamical instability is 
whether or not it is hysteretic. Our calculation has shown 
that the instability studied here is not, if a large magnetic 
field is applied. Without applied magnetic field, the na- 
ture of the spin wave instability depends on the precise 
form of the magnetic anisotropy, and both hysteretic and 
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non-hysteretic behavior can be expected, in principle. 

A noteworthy aspect of our calculation is that the spin- 
transfer torque is calculated self-consistently: the mag- 
nitude and direction of the spin-transfer torque depends 
on the spin accumulation in the normal metal, which, 
in turn, depends on the precise magnetization profile of 
the ferromagnet. In doing this, our work connects the 
the circuit theory for hybrid ferromagnet-normal-metal 
systems, which has been used extensively to describe 
the magnet's effect on spin accumulations in macrospin 
approximation^ and micromagnetic simulations, which, 
to date, have been restricted to simplified models for the 
spin-transfer torque. However, our simulations should be 
considered a proof-of-principle. They lack the spatial res- 
olution and sophistication that full-scale micromagnetic 
simulations have. 
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APPENDIX 

The perturbative calculation of Sec. focused on the 
case of a large applied magnetic field because, in that 
case, theoretical results for the spin wave amplitude do 
not depend on sample-dependent anisotropy energies. In 
this appendix wc outline the theory for the general case. 

For the most general case, one needs a better ansatz 
for the intrinsic torque Tan than Eqs. (|14|l or (|15|l . as well 



as an expression for the current-induced torque that does 
not rely on rotation symmetry around the easy axis. The 
general expression for the torque Tan is most conveniently 
derived from the Free energy, Tan = ^ (7 /M) (dF/dm) x 
m. Since we are interested in the mode m(qc) only, we 
can expand F in powers of mi(qc) and m2(qc), up to 
fourth oder as 

F{m) ^ ^ [A:imi(qc)^ + fc2TO2(qc)^] 

3 

-F^fc^mi(qc)^m2(qc)^"^' 

3=0 
4 

+ ^fc>i(qc)^m2(qc)^-^'. (55) 

3=0 

The higher-order expansion constants k'j, j = 0,1,2,3, 
and fc", j = 0,1,..., 4, are not governed by any spe- 
cial symmetry and therefore likely to be sample specific. 
(The cubic terms in the expansion of F{ui) may be for- 
bidden if there is a reflection symmetry around the easy 
axis. However, there is no such symmetry in the pres- 
ence of an applied magnetic field that is not aligned with 
the one of the sample's easy or hard axes, so that cubic 
terms need to be included in a general treatment.) Note 
that the higher-order torque terms are as important in 
determining the spin wave amplitude as the higher-order 
current-induced spin-transfer torque. Unless these coeffi- 
cients are known independently, a calculation of the spin 
wave amplitude has no predictive value — that was the 
reason why the main text addressed the case of a large 
applied magnetic field. 

We now list our general results for the second and third 
order potentials and third-order spin-transfer torque. 
The symbols used are defined in Sec.^of the main text. 
The second-order charge potential expansion coefficients 
for the normal-metal spacer are 



(2). ^ _ .y^ f [Gs(0) + G.(q) - 2Gs{q')]D{ci') + 2[Gs(0) - a(q')][Gs(q) - Gs(q')][gi + Gs(q')] 
q' 

X [TOi(q')r7ii(q - q') + TO2(q')TO2(q - q')] 

2g2[Gs(0)-a(q)][Gs(q)-a(qO] , , , , „ , r,,\ 
P^-Tn TTTn 7-\ [mi(q)m2(q-q) - m2(q)mi q-q) \ 

^(q ).9m(0)gm(q) J 

s / .92[Gs(q) - Gs(q')]Gs(q) , , i\ ■ , i\ , , i\ ■ 1 |\^ 

-^2^ I i:>(qOgm(q) - q )™i(q ) + ™2(q - q )™2(q )] 

q' " 

, [G.(q) - G.(qO][J(qO - G.(qO(ffi -|- G.(qO)] ^ r, ■ , r,\ 

+ \mi (q-q )m2 (q ) - ms (q - q )mi (q )] | . (56) 

f 2) 

The coefficient Oc (0) determines how the spin wave instability affects the resistance of the ferromagnetic layer, cf. 
Eq. 144|) in the main text. The second order correction to the 3-component of the spin potential is given by the 
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expansion coefficients 

as3^(q) = ej^|^^^^[mi(q')mi(q-q') + m2(q')TO2(q-q')] 
q' ™ 

- n( n L ( ^ + 2Ge(q))(Gs(0) + G.(q) - 2G.(q')) 

-D(q')3m(0)5m(q) 

- (G,(0) - G,(q'))(gi + G,(q))(.gn.(q) - 2{g+ + 2Gc(q)))(G,(q) - G,(q'))] 
X [TOi(q')TOi(q - q') + m2(q')TO2(q - q')] 

g2[GM - Gs(q)][2(g+ + 2Gc(q))(G.(q) - G.(qO) ^ gM] r , ^ . . ^ . 

Fv-7^ 7?v^ — TTn mi(q 7712 q~q -m2 qjmi q-q ^ 

^(q ).9m(0)gm(q) J 

- ( ^-^-(-^')[-'^-(-^) - + 2G (q))(G.(q) - G.(qO)] ^^^^^ _ ^ ^^^^ _ 

, I 2D(q'),gm(q) 
q' 

[Diq') - G,(q')(5i + Gs(q'))] [5m(q) - 2{g+ + 2Ge(q))(Gs(q) - G,(q'))] 



2i?(q')5m(q) 

ii(q - q')"i2(q') - m-2(q - q')mi(q')]| , 



(57) 



The very first term describes the effect of a uniform magnetization rotation; the remaining terms are the result of 
a non-uniform magnetization. There are second-order corrections to the spin potential expansion coefficients ai and 
02 that arise from the presence of cubic terms in the anisotropy Free energy. Such cubic terms cause second-order 
contributions to the time derivatives rhi and m2, which give a contribution to the second order spin potentials a^^-* 
in the same way as the first-order time contribution to the time derivative affects the first-order spin potentials a^^\ 
see Eq. (jHl- 

Instead of listing the third-order potentials a*^^^ , we describe the corresponding current-induced torque. We specialize 
to the contributions that are cubic in the magnetization amplitude at wavevector qc. The resulting torque has terms 
proportional to the third-order contributions to the time derivatives of the magnetization. These terms give rise to a 
renormalized Gilbert damping coefficient and a renormalized gyromagnetic ratio, see Eq. 122|) . The remaining terms 
can be written as 2r(0) -I- T(2qc), where (again using two-component spinor notation) 

f(3)(k) "^^ 



2Mde^D{q,)^g^{0)gM 
X {2ejg2i5(qc)Gs(qc)[Gs(qc) - Gs(k)][.g+ + 2Gc(k)][Gs(0) - G,(qc)]ia2m(qc)mT(qc)TO(qc) 
- [5,„(k) + 2(g+ + 2Gc(k))(G,(qc) - Gs(k))] 

X [ - ej[(2G,(qc) - Gs(0) - Gs(k))i?(qc) + 2(G,(0) - Gs(qc))(.9i + Gs(qc))(Gs(qc) - G,(k))] 

X [D{qc) - Gs(qc)(5i + Gs{qc))]m{qc)m^ {q.c)m{qc) 
+ ej52Gs(qc)[Gs(qc) - Gs(k)] [i?(qc) + 2(Gs(0) - Gs(qc))(5i + Gs{qcW<J2m{qc)m^ {q,)m{qc) 
+ %2Gs(qc)gm(0)[Gs(qc) - Gs(k)][D(qc) - Gs(qc)(.9i + Gs(qc))]TO(qc)TO'^(qc)"i(qc) 
+ (;i/2)[Gs(qc) - Gs{k)]Diq,)g^iO)[D{q,) - G,(qe)(2gi + Gs(qc))] 

X [m(qc)m'^(qc)«cr2"i(qc) - «cr2'7i(qc)m'^(qc)m(qc)] 
+ (V2)[Ge(qc) - Gs(k)]g,„(0)[(I?(qc) - Gs(qc)(ffi + G,(qe)))' - {g2Gsiqc)f]^a2m{qc)m^i^c)m{q,) ] } . 

(58) 



Once the perturbative expansions for the anisotropy known, the Landau-Lifshitz-Gilbert equation can be 
torque and the current-induced spin-transfer torque are solved for the magnetization dynamics. 
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